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摘要 采用 6250 型 近 红外 光谱 分 析 仪 ， 从 波长 680 到 1235 nm 分 别 对 日 系 乔 肾 217 颗 、 中 系 
du 174 颗 、 杂 交 种 看 肾 330 颗 做 了 扫描 测试 ， 结 果 表 明 : 上 肉 雄 乔 肾 的 近 红 外 光谱 有 极 显著 
的 差异 。 从 日 系 217 颗 乍 晴 中 随机 选取 雌雄 各 50 颗 作 为 学 习 样 本 ， 采 用 逐步 判别 方法 从 其 二 
阶 导数 光谱 中 抽取 5 个 特征 波长 作为 观测 向 量 ， 利 用 二 次 型 判别 函数 对 其 余 621 颗 蛋 里 的 性 
别 进行 识别 ， 其 准确 率 达 98.7%。 该 方法 明显 优 于 以 看 肾 重 量 或 大 小 识别 看 蜂 肉 雄 的 方法 。 
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家 在 ( Bombyx mori) 肾 的 肉 雄 上 自动 识别 研究 是 答 业 工作 者 长 期 想 解决 的 课题 之 一 ， 
因为 它 会 给 家 看 良种 繁育 和 育种 带 来 极 大 方便 ， 具 有 重要 的 经 济 意 义 。 家 看 遗传 学 家 利 
用 染色 体 易 位 ， 育 成 了 一 些 限 性 品种 趾 ， 这 对 区 分 看 或 章 的 肉 雄 是 一 个 突破 性 进展 ， 但 
由 于 染色 体 易 位 的 影响 ， 限 性 品种 在 某 些 经 济 性 状 方面 ， 还 赶不上 非 限 性 品种 ， 所 以 生 
产 上 目前 主要 以 推广 非 限 性 品种 为 主 。 刘 敬 权 等 中 利用 紫外 线 激发 严 光 识别 看 草 上 肉 雄 ， 但 
这 种 方法 同样 受到 所 谓 “ 判 性 品种 ”的 限制 ， 并 且 草 过 奖 光 的 遗传 规律 尚未 搞 清 楚 。 

近年 来 ， 近 红外 光谱 分 析 技 术 已 广泛 应 用 于 工业 ,农业 ,医药 等 各 个 方面 S44， 这 是 
一 种 非 破坏 性 的 仪器 分 析 法 ， 具 有 安全 可 靠 ， 无 致 变 作用 的 特点 ， 对 家 看 育种 和 良种 繁 
育 中 的 雌雄 鉴别 比较 合 通 。 本 研究 首次 将 近 红 外 光谱 分 析 技 术 和 模式 识别 理论 相 结合 ， 
将 其 应 用 于 和 寿 皇 和 晴 的 模式 识别 5， 试 图 为 鼻音 和 是 的 雌雄 识别 目 动 化 ， 从 理论 上 和 方 
法 上 打下 基础 。 这 里 仅 对 上 肉 雄 看 肾 近 红外 光谱 的 差别 及 其 识别 方法 作 一 详细 报道 。 


1 材料 和 方法 
1.1 实验 材料 及 测试 
本 试验 所 用 的 乔 肾 ， 夏 乔 期 为 杂交 种 ， 由 徐州 市 草 丝 绸 研究 所 提供 : 中 秋 和 看 期 为 中 
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系 品种 和 日 系 品 种 ， 由 江苏 省 镇 江 答 种 场 选 种 站 提供 。 测 试 工 作 在 北京 市 蔬菜 研究 中 心 
生化 分 析 室 (农业 部 开放 实验 室 ) 进行。 测试 仪器 为 美国 太平 详 公 司 生产 的 6250 型 近 红 
外 光谱 分 析 仪 和 D 一 744 型 光 导 纤维 探头 ， 光 栅 式 扫描 ， 波 长 范围 680 ~ 1235nm, Hi 
硅 检测 器 测量 反射 信息 。 被 测 乔 里 发 育 进 度 从 复眼 着 色 到 触角 着 色 之 间 ，。 
首先 前 草 取 肾 ， 对 每 颗 看 肾 做 近 红 外 扫描 ， 由 计算 机 自动 将 反射 信息 转换 成 
log (1/4) 值 贮存 ， 其 次 称 肾 重 ， 鉴 别 下 肾 肉 雄 。 本 试验 于 1994 FE E BH UM A RACE 
44 330 颗 ， 其 中 肉 肾 150 Ei, MENS 180 Wi; 于 1994 年 中 秋 期 测试 中 系 乔 肾 174 S, H 
"Htm 86 颗 ， 梭 肾 88 颗 ， 测 试 日 系 乔 肾 217 颗 ， 其 中 肉 肾 100 88, 28598 117 颗 。 


1.2 ” 肉 雄 每 晴 的 识别 方法 


采用 NSAS 光谱 数学 程序 ， 求 原始 光谱 (log 1 / R) 的 二 阶 导 数 光 谱 。 由 学 习 样 本 二 
阶 导 数 光 谱 数 据 ， 采 用 逐步 判别 方法 ， 选 取 与 重 肾 性 别 有 关 的 特征 波长 上 。 比 较 肉 雄 看 里 
光谱 特征 的 差别 ， 并 对 其 做 统计 检验 。 由 学 习 样 本 的 人 选 特征 建立 值 Bayes 二 次 型 判别 
PREX: 

u, Q)-(x-x"ysx-x?)«mnlSI-2ing i=l, 2 

AP, X= xs 5s x. x” j=l, 2.5 p A59 SR BMEHETE XI XR p 个 人 选 特 
Efl, X 0, S. q OA ARRERA RERE, A EA ERAR. MERDA 
B XPEZE—AGKAUEZAK XAORA HAREMA: Eu OO «u, OO WA AE, RA 
X HE. 

判别 效果 的 检验 方法 : (DIDE: 从 于 个 学 习 样 本 中 逐一 剔除 一 个 样本 ， 以 n 一 1 
个 学 习 样 本 建立 判别 函数 ， 对 剔除 样本 做 判别 。 求 对 关 个 噜 除 样本 的 错 判 率 .(2) 用 学 习 
样本 建立 判别 函数 对 检验 样本 预测 。 以 上 判别 分 析 均 在 Quick Basic 环境 下 编程 ， 计 算 
程序 设计 采用 双 精 度 变量 。 


2 结果 与 分 析 
2.1 ”次 奴 大 晴 近 红外 扫描 光谱 的 差异 及 统计 特征 


1.3 - D b 
E x 0.10 9 
= oq = 
È 2 | 
a i i á 
入 Æ —0.10 
Eu 让 
Ws 0.9 





900 1 000 1 100 1 200 900 1 000 i 100 1 200 
i& K (nm) 波长 (nm) 


图 1 MERA 50 个 样本 的 半 均 扫描 光谱 
a， 示 原始 光谱 log (1 /R); b， 示 二 阶 导数 光谱 
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M.F AR ER 217 个 样本 中 随机 选取 肉 雄 各 50 个 ， 分 别 求 其 原始 光谱 和 二 阶 导数 光 
谱 的 平均 数 ， 由 于 雌雄 蚕 晴 的 光谱 差异 主要 反映 在 900 ~ 1 200 nm 的 波段 上 ， 图 1 仅 
绘制 了 这 个 波段 上 的 光谱 图 。 

比较 平均 原始 光谱 和 平均 二 阶 导 数 光谱 ， 似 乎 原始 光谱 更 能 反映 峻 雄 间 的 差异 。 其 
实 不 然 ， 因 为 看 肾 原始 光谱 的 差别 ， 受 乔 肾 重量 (大 小 ) 的 影响 较 大 ， 和 看 肾 越 大 ， 反 射 的 
能 量 也 就 越 多 ， 但 是 鼻 肾 重量 在 个 体 间 的 变异 较 大 ， 从 后 面 的 分 析 可 以 看 到 ， 仅 依 乔 肾 
重量 区 分 看 肾 肉 雄 ， 错 判 率 是 很 高 的 。 一 ,二 阶 导数 光谱 反映 的 是 原始 光谱 的 变化 率 ， 
它们 的 差别 与 晴 体 的 理化 性 质 有 关 ， 虽 然 在 图 示 的 精度 上 ， 反 映 的 光谱 差异 较 小 ， 但 
是 ， 二 阶 导 数 光 谱 变 异 较 小 ， 比 较 稳定 ， 更 能 反映 出 肉 雄 间 的 差异 。 

由 于 在 扫描 区 间 上 的 光谱 值 达 700 个 ， 为 简化 分 析 ， 采 用 逐步 判别 方法 ， 从 众多 的 
波长 测试 数据 中 抽取 与 蛋 肾 性 别 关系 密切 的 特征 波长 5 个 ， 现 将 5 个 波长 处 的 原始 光谱 
值 和 二 阶 导数 光谱 值 分 别 列 于 表 1 ME 2. 


表 1 原始 光谱 在 5 个 特征 波长 处 的 统计 特征 


波 长 (nm) 921 933 1041 1 099 1182 
BE d 平均 数 0.94076 0.94401 0.95996 0.89324 1.2448] 
标准 差 0.04277 0.04244 0.04001 0.03992 0.03797 

Ht own 平均 数 1.03666 1.04017 1.04793 0.98358 1.33610 
标准 差 0.04091 0.04055 0.30846 0.03830 0.03850 

F- fü 102.29" 115.307 92.36" 13.78" 8.71" 


注 : 下 检验 为 附加 信息 法 ，” 


表示 1% 的 显著 水 平 ， 瑟 ,(1，94) =6.91， 下 同 
x2 二 阶 导数 光谱 在 个 特征 波长 处 的 统计 特征 


波 K (nm) 921 933 1041 1 099 1182 
HE dg 平均 数 0.02435 0.04861 0.01229 0.02407 —0.15727 
标准 差 0.00115 0.00238 0.00097 0.00100 0.00363 

M 肾 平均 数 0.02251 0.04281 0.01163 0.02176 — 0.16301 
标准 差 0.00113 0.00215 0.00069 0.00099 0.00344 

F- fü 195.70" 319.29" 58.71" 43.41" 29.29" 


表 中 5 T CI ABS BCIBOGTA HU ZEB CERERI ARORURCBUM 0.937, 55 SE 88 PE RECO 和 
1) 的 复 相 关系 数 为 0.925， 而 二 阶 导 数 光谱 值 与 晴 重 的 复 相关 系数 为 0.935， 与 蛋 肾 性 别 
的 复 相关 系数 为 0.972， 可 见 ， 二 阶 导 数 光谱 更 能 反映 肉 雄 看 肾 间 的 差异 。 


2.2 ”上 典 雄 乔 肾 的 光谱 识别 
以 上 面 日 系 品 种 的 100 个 样本 的 测试 数据 作为 学 习 样 本 ，921，933，1 041，1 099， 
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1 182 nm 5 个 波长 处 的 二 阶 导数 值 作为 观察 向 量 ， 建 立 Bayes KWAP RA, KAJ 
切 法 估算 错误 率 ， 并 以 此 判别 函数 对 其 余 117 个 日 系 品 种 的 光谱 进行 识别 ， 其 识别 效果 
见 表 3. 

在 100 个 学 习 样 本 的 刀 切 法 检验 和 对 117 个 检验 样本 的 预测 中 ， 仅 有 2 个 错 判 。 根 
据 试 验 计 录 知 道 其 中 一 个 是 僵 死 是 ( 雌 ) ， 另 一 个 是 软化 的 不 正常 晴 ( 雄 ) 。 为 进一步 检验 
对 不 同 品种 雌雄 晴 的 识别 效果 ， 我 们 仍 用 前 面 所 建立 的 二 次 型 判别 函数 ， 对 174 个 中 系 
TRI 330 个 杂交 种 看 肾 进行 肉 雄 识别 ， 结 果 见 表 4。 


表 3 二 次 型 判别 函数 对 日 系 每 旺 表 4 二 次 型 判别 函数 对 不 同 品 种 
的 雌雄 识别 效果 簿 师 的 雁 雄 识别 效果 
判别 结果 判别 结果 
OO l| qw e| 一 一 一 一 一 一 一 
RUM 方法 1 方法 2 方法 1 方法 2 
HO Hu hi HO i H ig 
i 30 0 49 ] hn 86 0 149 1 
ME 0 50 | 66 雄 0 88 5 175 
错 判 率 ( %) 0 1.71 错 判 率 ( %) 0 1.82 
方法 k 100 个 学 习 样本 的 刀 切 法 检验 ; 方法 1， 对 中 系 品种 174 个 检验 样本 的 预测 : 
方法 2 对 日 系 品种 117 个 检验 样本 的 预测 方法 2 对 杂交 种 330 个 检验 样本 的 预测 


虽然 中 系 品 种 .日 系 品种 和 杂交 种 的 蚕 晴 在 形态 和 重量 上 有 一 定 的 差别 ， 测 试 时 间 
也 不 尽 相 同 ， 但 品种 间 的 上 肉 雄 识别 仍然 取得 了 很 好 的 效果 ， 对 621 个 检验 看 肾 性 别 预测 
的 准确 率 达 98.7%， 这 也 说 明 不 同 品种 的 二 阶 导数 光谱 在 反映 峻 雄 晴 的 差异 方面 具有 共 
同 特征 。 
2.3 ”以 晴 重 为 特征 值 的 肉 雄 判别 效果 与 光谱 识别 的 比较 

人 们 曾 以 牌 昭 重 量 来 区 分 雌雄 短 晴 ， 其 效果 如 何 ? 为 与 光谱 识别 法 作 一 比较 ， 我 们 
仍 用 100 个 日 系 鼻 晴 为 学 习 样 本 ， 以 等 晴 重 量 为 特征 值 建立 判别 函数 ， 对 其 余 621 个 日 
系 、 中 系 、 杂交 种 鼻 肾 做 上 肉 雄 鉴别 ， 其 结果 见 表 5。 


比较 表 3, 3€ 4、 表 5 可 知 ， 依 肾 重 识别 表 5 以 肾 重 为 特征 值 对 系 师 的 肉 雄 判别 效果 
看 肾 上 肉 雄 的 错 判 率 是 光谱 识别 法 的 12 fü. = 一 3 


这 是 由 于 看 肾 重 量 在 品种 内 和 品种 间 有 较 — DUE 

大 变异 ， 造 成 肉 雄 肾 间 的 差异 不 明显 之 缘 方法 1 方法 2 

故 。 我 们 还 知道 ， 乔 肾 重 量 与 蚕 期 的 饲养 E E E g 
环境 条 件 有 很 大 关系 ， 即 使 是 同一 品种 ， H 6 7 22 D 
在 不 同 环境 下 饲养 ， 肾 重 也 会 有 很 大 差 雄 2 48 76 259 
55. BREL AK SR RECRCRE JE DX A RE BER SE 错 判 率 ( %) 9 15.46 


NUBE TE, 而 用 和 蛋 师 的 光谱 识别 。 方法 1 100 个 学 习 样本 的 刀 切 法 检验 
方法 ， 基 本 上 可 以 排除 由 于 品种 的 不 同和 HEL 对 日 系 .中 系 .杂交 种 621 个 检验 样本 的 预测 
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3 Wit 


影响 雌雄 重 晴 近 红 外 光谱 差异 的 原因 是 值得 进一步 研究 的 。 这 种 差别 并 不 是 某 些 品 
种 特有 的 现象 ， 因 为 我 们 选用 的 中 系 . 日 系 以 及 杂交 品种 都 是 从 生产 上 随机 选取 的 ， 它 
们 的 二 阶 导数 光谱 都 表现 了 共同 的 特征 。 如 果 说 这 种 差别 来 自 于 肾 体 的 大 小 或 重量 ， 从 表 
5 已 清楚 地 看 到 ， 仪 由 肾 重 是 难以 识别 乔 肾 肉 雄 的 。 我 们 从 入 选 的 5 个 二 阶 导 数 特征 值 
与 肾 重 的 相关 系数 (0.935)， 和 与 看 肾 肉 雄 的 相关 系数 ( 0.972) 来 看 ， 入 选 的 特征 值 中 除了 
包含 肾 重 的 信息 外 ， 还 一 定 包含 其 他 与 肉 雄 性 别 有 关 的 信息 ， 我 们 初步 推测 这 与 肾 体 内 
看 卵 的 有 无 ， 脂 肪 的 含量 以 及 含水 率 有 一 定 的 关系 ， 虽 然 我 们 还 没有 进一步 证 实 ， 但 从 
近 红 外 技术 的 其 它 应 用 领域 来 看 * 4， 我 们 可 以 概括 地 说 ， 答 晴 近 红外 光谱 的 差别 是 是 
体 的 大 小 和 体内 化 学 成 分 的 综合 影响 ， 

表 1 和 表 2 中， 虽然 人 选 特征 值 的 平均 数 在 肉 雄 群体 间 存 在 极 显 着 的 差异 ， 但 从 每 
个 特征 值 的 标准 差 来 看 ， 仅 仅 依 靠 某 一 特征 值 是 难以 区 分 看 肾 肉 雄 的 ， 也 就 是 说 肉 雄 乔 肾 
的 光谱 差异 是 光谱 整体 结构 的 综合 差异 ， 而 不 是 某 一 二 个 波长 特征 所 能 描述 的 。 乔 肾 
光谱 识别 的 任务 ， 就 是 通过 适当 的 数学 手段 ， 使 其 光谱 的 结构 能 够 反映 出 雌雄 间 的 差 
异 ， 并 从 众多 的 光谱 数据 中 ， 抽 取 少 数 与 看 肾 性 别 有 关 的 特征 值 ， 根 据 这 些 特征 值 构造 
判别 函数 ， 综 合 地 反映 肉 雄 看 肾 光 谱 结 构 的 差异 ， 从 而 达到 上 峻 雄 识 别 的 目的 。 因 此 ， 光 
谱 识 别 的 效果 关键 在 于 我 们 能 否 利用 光谱 中 所 包含 的 有 用 信息 。 许 多 看 业 工 作者 也 曾 利 
用 射线 或 光谱 识别 看 肾 或 看 草 的 肉 雄 ， 其 未 成 功 的 原因 很 可 能 就 是 没有 对 数据 进行 合理 
的 分 析 和 利用 ， 因 为 光谱 的 微小 差别 有 时 只 赁 肉眼 难以 察 党 ， 通 过 适当 的 数学 变换 ， 就 
可 以 把 这 种 差别 反映 出 来 。 本 人 研究 正 是 结合 了 模式 识别 的 原理 和 近 红 外 光谱 分 析 技 术 ， 
才 达 到 了 对 看 肾 光 谱 识 别 的 效果 。 

当然 ， 本 研究 只 是 对 蛋 肾 光谱 识别 方法 和 可 行 性 的 探讨 ， 离 蛋 肾 的 自动 识别 还 有 一 
段 距离 ， 还 有 许多 问题 值得 进一步 研究 。 辟 如 学 习 样 本 的 代表 性 和 特征 值 选 取 的 合理 性 
问题 ， 可 以 预计 ， 当 这 些 问题 解决 之 后 ， 和 看 肾 识 别 的 准确 率 还 能 进一步 提高 。 当 我 们 通 
过 大 量 的 实验 确定 了 少数 几 个 波长 之 后 ， 对 看 肾 的 光谱 识别 ， 就 不 需要 做 680 ~ 1 235nm 
的 全 波段 扫描 了 ， 这 样 可 以 大 大 减少 测试 时 间 ， 
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DIFFERENCE BETWEEN MALE AND FEMALE PUPAE OF 
BOMBYX MORI IN NEAR INFRARED REFLECTANCE 
SPECTRA AND THEIR PATTERN RECOGNITION 
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Uniwrsity Xuzhou 221009) Zhenjiang 212018) 
Sun Aiqun Jin Tongming 
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institute Xuzhou 221009) Beijing 100081) 


Abstract In this paper, difference between male and female pupae of Bombyx mori 
in near infrared reflectance spectral NIRS) and the feasibility and method of spectra 
recognition. to identify the sex of fresh pupae were studied. By Model 6250 NIR 
instrument, indestructive determination of 721 fresh silkworm pupae selected from 
Chinese race, Japanese race and their hybrid was made with suitable wavelengths 
chosen between 680nm and 1235nm. Fifty male and fifty female pupae was randomly 
selected from 217 pupae of Japanese race as learning samples. The feature at 5 
wavelengths was extracted as observation vector from data of their second derivative 
spectra by a stepwise discriminant analysis. Statistical analysis showed that there were 
extremely significant differences between male and female pupae In NIRS at the 5 
wavelengths. When the Bayes quadratic discriminant function built by the feature was 
used to predict the sex of the other 621 test samples, 98.794 of correct rate of the 
identification recognition was obtained. This method is much better than other 
identification ones which use weight or size of pupae as their feature. The mechanism 
is discussed. 


Key words silkworm (Bombyx mori) pupa, determination of pupa sex, NIRS, pattern 
recognition. 


